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I. INTRODUCTION 

Vulcanization has a vast array of technological 
and commercial applications. The vulcanization pro- 
cess transforms a macromolecular liquid into an amor- 
phous solid by randomly introducing permanent (chem- 
ical) crosslinks. In the liquid state all macromolecules 
are delocalized. In the amorphous solid state a non-zero 
fraction of macromolecules forms a macroscopic network. 
The constituents of this network are localized at random 
positions about which they execute thermal motion char- 
acterized by a distribution of finite localization lengths. 
The vulcanization transition (VT) between the two states 
(occurring at a critical density of crosslinks) represents a 
continuous phase transition. 

Due to the work of Goldbart, Zippelius, and coworkers 
(GZ) 1^ a rather comprehensive theoretical description of 
the VT exists to date on the level of a mean field approx- 
imation. The mean field theory gave a first glance on the 
relation between the VT and the percolation transition, 
and the critical exponents describing the VT were shown 
to be consistent with the mean field critical exponents for 
percolation. Peng et al. Q introduced a minimal model 
for the VT and discussed it under the aspects of uni- 
versality as a common theoretical formulation of general 
amorphous solidification transitions (of which the VT is 
a prime example). Recently Peng and Goldbart (PG) 
carried out a renormalization group improved one-loop 
calculation based on the minimal model |^ . Their calcu- 
lation showed that the critical exponents of vulcanization 
and percolation are in conformity to one-loop level. 

In the past, the VT has often been addressed directly 
from the perspective of percolation theories |^ . In con- 
trast to the work of GZ this perspective takes into ac- 
count only a single ensemble of random connections, 
and therefore does not incorporate the effects of both 
quenched randomness and thermal fluctuations. Given 
that an essential aspect of the VT is the impact of the 
quenched random constraints on the thermal motion of 
the constituents, the a priori identification of the VT 
with percolation is a nontrivial matter and one may ask: 



What is the relevance of percolation theory to the VT? 

The purpose of this paper here is to explore the connec- 
tion between vulcanization and percolation in depth. We 
compare the minimal model |^ to a field theoretic model 
for the random resistor network (RRN) which we recently 
studied . Such a RRN is nothing more than a percola- 
tion model where randomly occupied bonds are assigned 
a finite nonzero conductance. In contrast to PG, which 
implement a "momentum shell" renormalization group, 
we use the powerful methods of renormalized local field 
theory. The momentum shell procedure, although phys- 
ically intuitive, mixes scaling properties with features of 
other universal quantities like scaling functions. There- 
fore it is extremely difficult to extend the procedure be- 
yond first order in perturbation theory. Renormalized 
field theory goes the other direction. Instead of the elim- 
ination of fluctuating degrees of freedom near the upper 
momentum cutoff and a rescaling of the momenta, the 
momentum cutoff is send to infinity at the beginning. 
The ultraviolet (UV) infinities resulting from this limit- 
ing procedure are eliminated by the so called renormal- 
ization factors. These renormalization factors then deter- 
mine the Gell-Mann-Low renormalization group equa- 
tion (RGE) which encodes all the scaling properties of 
the theory. This RGE represents the cleanest way to an- 
alyze the full scaling structure and its critical exponents. 
Physical properties beyond the scaling structure, how- 
ever, cannot be inferred from the RGE. Indeed, these 
may be distinct for different theories like vulcanization 
and percolation. A further advantage of renormalized 
field theory is that it allows to find out general proper- 
ties of the renormalization factors (and hence the RGE) 
to all orders of perturbation theory because they result 
solely from so called superficially divergent Feynman di- 
agrams. Our analysis presented in this paper thrives on 
this advantage. We show rigorously that vulcanization 
and percolation involve the same primitive divergences 
and hence are governed by the same critical exponents. 
Furthermore, we compare the order parameters of the 
VT and the RRN. 
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II. MODELLING VULCANIZATION AND 
PERCOLATION 

The theory of GZ and PG is based on the order pa- 
rameter field 



ik-f 



(1) 



defined on the replicated d-dimensional real space with 
coordinates f = (ri,--- ,r„). The corresponding repli- 
cated wave vectors are k — (ki, • • • ,k„). Note that our 
n corresponds to (n-l- 1) in the theory by GZ. The replica 
limit n ^ 1 has to be taken before any other limit. The 
volume V of the real space is considered as being finite. 
Thus, fc is a discrete vector. • • • is an abbreviation for 
• • • which is in the infinite volume limit equivalent 
to (27r)~"'' J d"''^k ■ ■ ■ . Microscopically, the order param- 
eter is a n-fold correlation function of density fluctua- 
tions and characterizes the amorphous state. The fluc- 
tuations about the average density thcm-self represent a 
non-critical stochastic variable which is excluded by the 
constraints 



Jl d'^r^^'^ n(f) = , a = I,-- - ,n . (2) 



Hence, Vl{k) is only non-zero if k belongs to the "higher 
replica sector" (HRS), that means if at least two distinct 
wave vectors k("),k(^) of k are non-zero. The minimal 
model for the vulcanization transition O] is defined by 



WvT = / d 



Mr? 



9 



(3) 



where t—Tc measures the distance to the critical crosslink 
density at Tc- In a mean field approximation r is positive 
in the liquid phase and the VT occurs aX t — Tc = Q. 
By virtue of the usual identification of ^{f) with the 
microscopic density correlation, g is a positive coupling 
constant. The Hamiltonian TivT is complete in the renor- 
malization group sense, i.e., it contains all relevant cou- 
plings and neglects all irrelevant ones. 

We are going to compare the perturbation theory based 
on the minimal model TivT to the perturbation theory 
of the field theoretic model for the RRN §, 0, §. This 
field theoretic model is based on an order parameter field 
(/)(x, 0) which lives on the d-dimensional real space de- 
noted by the coordinates x. The variable 9 denotes the 
_D-fold replicated voltage at position x. For regulariza- 
tion purposes, 9 = vA9 takes discrete values on a D- 
dimensional torus, the replica space, i.e., v is chosen to 
be a /^-dimensional vector with integer components i^*-"-* 
satisfying -M < < M and z^(") = mod (2A/). 

The order parameter field is restricted by the condition 



(4) 



It follows that the replica space Fourier transform (^(x, A) 
of the order parameter field, defined by (/)(x, 9) = 
{2M)-° J2x "^(x, A) exp(iA • 9), satisfies <^(x, A = 0) = 0. 
Without exception we study the limit _D — > 0, A/ ^ cx) of 
the rephca space with (2M)^ 1 and A9 = 9o/VM 
0. Here 6*0 is a constant which sets the width of the volt- 
age interval such that [-9oVM < 9^"^ < 9oVM]. In 
the limit D ^ 0, M — > oo the constant 6*0 plays the 
role of a redundant scaling parameter, i.e., the theory 
is independent of its value. In the following we write 

Er ■ ■ (A^)^ Ee---^ Wd---=;- le---. where the 
approximations become exact in the limit studied. 
The field theoretic Hamiltonian for the RRN reads 



Hrrn = / d'^x 



^0(x,0y^i [v0(x,r'' 



w 
2" 
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x,( 



(5) 



The parameter r — ^ {pc —p) specifies the deviation of 
the occupation probability p from the critical probability 
Pc- In mean field theory the percolation transition hap- 
pens at T = Tc = 0. w is proportional to the resistance 
of the individual random bonds. 

In the percolating phase, t < Tc, the mean order pa- 
rameter is given by 



(x,A) 



\2 

cxp( - Yi?oo(x) 



C 



Xoo(x) exp ( - —R^{yi) 



Poo ( exp 



(6) 



Here i?oo(x) is the (random) resistance between an arbi- 
trary point X and infinity. Xoo(x) is an indicator func- 
tion which is one if x is connected to infinity, i.e., if x 
belongs to a percolating infinite cluster, and zero other- 
wise. (• • • )c denotes the disorder average over all config- 
urations of the diluted lattice. (• • • )p stands for disorder 
averaging conditional to the constraint that x belongs to 
an infinite cluster. Note that Poo = (Xoo(x))c is the per- 
colation probability, i.e., the order parameter for usual 
(purely geometric) percolation. 

'^RRN reduces for — > to the Hamiltonian for 
the n — (2M)^-state Potts model with n 1 for 
D ^ 0. This is important because the Potts model 
is known to describe percolation in this limit The 
connection between TYrrn and the Potts Hamiltonian 
becomes evident by relabeling the n fields (f) by an in- 
dex a — 1, ■ ■ ■ ,n: 0(x, 0) — > <j>a{x.). The constraint (^) 
then reads J^a^ai^) — 0. Taking the usual base e-"\ 



i = 1 



- 1, with 



0, E„ 



(a) 



5a f} — ^/n, and upon defining (n — 1) in- 



dependent fields s by 0a (x) = Ei Si(x), we get the 
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Potts-Hamiltoiiian in the form 



^.,(x)2 + i[V..(x)]^ 



I Aijfcs,(x)sj(x)sfe(x)| , 



(7) 



where A, 



'^ijk = J2a usual "Potts tensor". 

Note that the introduction of the finite bond resistance 
w reduces the Potts symmetry, i.e., the full permuta- 
tion symmetry of the n fields s belonging to the basic rep- 
resentation of the symmetry group Sn, to translation and 
rotation symmetry in the D-dimensional replica space. 



III. RENORMALIZATION AND SCALING 

We proceed with a renormalization group analysis of 
the VT by employing standard methods of renormalized 
field theory |1C[| . The principal diagrammatic elements 
are easily gathered from TivT- First, we have the three- 
leg vertex g. Because the corresponding interaction is 
cubic, c?c = 6 is found to be the upper critical dimension 
of the VT. Second, it follows from the quadratic part of 
TivT that the principal Gaussian propagator is given by 

n 

G^°"^(^) = G{k){l - ^£,k(")^(") + - ^)\6} - (8) 

a = l 

where G[k) — {t + k'^)^^ . e'") is a n-dimensional vector 
whose a-th component is one and all other components 
are zero. Accordingly, k*-"-* = k ■ e^"^ is the component of 
k containing replica a. The part of G^°^^{k) embraced 
by the curly brackets ensures the constraint k E HRS. 
Equation (||) can be interpreted so that G^"^'^ (A;) decom- 
poses into three generic parts one of which is proportional 
to (n — 1) and hence can be discarded from the onset. 
We call the remaining propagators G^™{k) = G{k) and 
G''"^(k'"'') — G(fc)(5j, ^(q) .(a) amorphous and liquid prop- 
agators, respectively, because the HRS fields entering 
G™(fc) are capable of diagnosing amorphous solidifica- 
tion in contrast to IRS fields constituting the G'"^(k^"'). 
Moreover, we refer to G'"^(k^"-') as having color a. Due 
to the propagator decomposition each principal diagram 
decomposes into a sum of amorphous diagrams consist- 
ing of amorphous and liquid propagators. Note that such 
an amorphous diagram may feature liquid propagators of 
just one or of several colors. 

Now we are in the position to address the renormal- 
ization of the VT Hamiltonian. For this task it is suf- 
ficient to consider those Feynman diagrams contribut- 
ing to the superficially divergent vertex functions for 
n —^ 1. Thus, we restrict ourselves to one-particle irre- 
ducible diagrams with two or three amputated external 
legs. In the following we denote generic diagrams of this 
type by D2{K) and D^dK}), respectively. K stands 



for an external wave vector and {K} is an abbreviation 
for {Ki, K2, K3 = —Ki — 7^2}- Once the decomposi- 
tion of the bold diagrams into amorphous diagrams has 
been accomplished (i.e., the constraint k € HRS is imple- 
mented) it is save to take the continuum limit in fc-space. 
Then we extract the divergences of these diagrams by ex- 
panding in a Taylor series in the external wave vectors, 
2)2(^) = ©2(6) +D^(6)i^2 and ^:i{{k}) = 2)3({6}), 
where higher order terms are discarded since they are 
superficially convergent. 

In the following we use without exception dimensional 
regularization. Thus, UV-divergencies appear as poles 
in the deviation e = 6 — c? from dc- These poles are 
eliminated from superficially divergent vertex functions 
by using the renormalization scheme 



{l = Z^^^n, T^T = Tc + Z-^ZrT , 
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(9) 



where Ag is a suitably chosen amplitude, is a con- 
venient length scale, and u is a dimensionless version of 
the coupling constant. In the minimal renormalization 
procedure, i.e., dimensional regularization in conjunction 
with minimal subtraction, fc is zero and the e-poles are 
eliminated by Z-factors of the form 



Z 



1 + E 



(10) 



m— 1 



The F'^'")(u) are expansions in the coupling constant u 
beginning with the power m™. A central theorem of 
renormalization theory, cf. Ref. jlOj, ensures that this 
procedure is suitable to eliminate the UV-divergencies 
from any vertex function order by order in perturbation 
theory. The Y^^^ (u) resulting from the primitive diver- 
gences determine the RGE entirely. 

First we analyze S)/({0}), I being any positive integer, 
under the assumption that it contains, apart from amor- 
phous propagators, liquid propagators of just a single 
color, say a. Obviously there are n different realizations 
of the generic type represented by Di{{0}) correspond- 
ing to the n different colors the liquid propagators can 
have. The sum over all these realizations appears in the 
decomposition of the pertaining bold diagram. However, 
all these realizations symbolize equivalent mathematical 
expressions since the internal wave vectors labeled by the 
colors are merely integration variables. Hence the overall 
contribution of this generic type to the primitive diver- 
gence of its bold diagram is 



(11) 



Next we assume that Di{{d}) contains liquid propaga- 
tors of several colors a, (3, - ■ ■ . Then the sum over the 
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realizations of this generic type is given by 



n n 



E E •••2)/({6}) 

= n(n-l)...S)i({0}), 



(12) 



i.e., it vanishes for n — > 1. 

Now we turn to 'D2{K) and assume that all featured 
liquid propagators are of the same color, a for sake. Since 
this color is distinguished, it is convenient to isolate it by 
setting k = i'a+e^^^k^") with = K-e(")k("^ Then 
D2{k) has the expansion D2{k) = £>2(0) + DqKI + 
^^jj-(")2 . . . _ It is important to realize that the expan- 
sion coefRcients D are independent of the choice of the 
color. Thus, the summation over all realizations of the 
generic type represented by D2{K) gives 



ES2W'^ 1)2(6) + 



(13) 



where we have exploited that J2a 



1 — 1)K^ van- 
ishes for n ^ 1. The case of several colors ••• 
can be analyzed as above by isolating the distinguished 
colors. D2{K) has the expansion '£)2{K) — ©2(0) + 



Dik(")2 + Dzk^^)^ + ■■■ with Kajj.^... 

K - e(")k(") - e^'^'k^''^ . Then it is straightforward 

to check that 



E E ■ 

a=l/3(#a) = l 



■^2iK) (x{n-l). 



(14) 



i.e., the contribution of the two- leg diagrams featuring 
liquid propagators of several colors vanishes for n ^ 1. 

To analyze the single color diagrams further, we 
write G{k) in Schwinger parameterization, G{k) = 
e^('^+'^^^*. Suppose that D2{k) comprises P prop- 
agators, being amorphous and P'"^ being liquid with 
color a. Assume that the diagram has L loops, i.e., we 
have to integrate over L independent combinations {qi}' 
of internal wave vectors {qi^ • • • , qp}. This integration is 
of the form 



/ expj- E UQh-T.'^^^"'"] 

•'{4.}' ifzp^rn jgp J 

X A({0)exp{-P({0)k(")2}, (15) 

where Qa.i ~ Qi — e'^'^^q'^"). {t}pam denotes the 
subset of the {ti^ ■ ■ ■ ,tp} belonging to the amorphous 
propagators. By comparing Eqs. ( p^ ) and ( p^ ) in 
the limit rt ^ 1 we see that 1)2(0) results from 
the integral J d{t}A{{t}) exp{—Tj2iti) and Di from 
J d{t}A{{t})B{{t})exp{-TY,.U). We conclude from 
Eq. ( p^ ) that these parts can be extracted directly by 



replacing all amorphous and liquid propagators by ele- 
mentary propagators of the type (r -I- q^)^^. Moreover, 
we can simplify K'^ to k'^ once the wave vector integra- 
tion has been carried out. An analogous, yet simpler, 
reasoning applies to DsdO}) containing a single color a. 
The case that D2{K) and SsdO}) contain solely amor- 
phous propagators, i.e., no colors at all, can be analyzed 
in a similar fashion. 

The quint-essence of our considerations is that prim- 
itive divergences steaming from diagrams with multiple 
colors drop out in the limit n — > 0. Hence it is suffi- 
cient for calculating the critical exponents of the VT to 
work with an arbitrary single color. Further, after the 
decomposition into amorphous diagrams all the propa- 
gators can be identified with elementary propagators of 
the type (T-|-q^)~^. We will see in the following that the 
perturbation series resulting from this effective decompo- 
sition coincides with the diagrammatic expansion of the 
field theory of RRN based on the Hamiltonian TYrrn in 
the limits D,w 0. 

The cubic interaction term in TYrrn leads to the three- 
leg vertex g. The principal propagator for the RRN is 
given by 

G'^°>'^(k,A) = G(k,A){l-5x,o}, (16) 

where G(k, A) = (t + k^ + wX'^y'^. Due to the fac- 
tor { 1 — (5^ which enforces the constraint A 7^ the 
principal propagator decomposes in a conducting part 
G'^°"'^(k, A) = G(k, A) carrying replica currents and an 
insulating part G'"^(k) ~ G(k, X)5^ ^ not carrying replica 
currents. Each principal diagram decomposes into a sum 
of conducting diagrams consisting of conducting and in- 
sulating propagators. As soon as this decomposition is 
accomplished it is save to switch to continuous replica 
currents. 

Now it is important to realize the one-to-one corre- 
spondence between conducting and amorphous propaga- 
tors as well as the one-to-one correspondence between 
insulating propagators and liquid propagators for n — > 1. 
As far as primitive divergencies are concerned, these 
one-to-one correspondences lead identical diagrammatic 
expansions (including the combinatorial fore-factors of 
the diagrams) for the VT and the RRN up to appar- 
ent distinctions in the propagators. Due to these dis- 
tinctions a diagram X)2(k, A) analogous to D2{K) (with 
pcond ^ pam) involvcs instead of Eq. (^ the integra- 
tions 

/ exp<^ E - E r 

= [F{{wt}p.a.a)]° exp{-C{{t}p.a.a)wX^} 
X A{{t})exp{-B{{t})k'}, (17) 
with the same functions A, B, C, F as in Eq. (p^). Equa- 



tion (|1J) leads for P) ^ to 

1)2 (k. A) 1)2(0) + DqwS? + Di-k^ + ■■■ , (18) 
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where the divergent coefficients 2)2(0) and Di are iden- 
tical to those appearing in Eq. (|3|). Here in the RRN, 
however, a further divergent coefficient Dq arises. Thus, 
the renormahzation scheme 



w ^ w — Z~^ZujW , r ^ f = f c + Z^^ ZtT (19) 

involves a further renormahzation factor Z^,. The other 
factors, Z, are identical to those in the renor- 

mahzation scheme (||) because the superficially divergent 
parts of the two diagrammatic expansions coincide to any 
order for w — 0. This leads to the conclusion that the 
renormahzation factors appearing in the renormahzation 
scheme for the VT are identical to those for percolation. 

Due to the independence of the unrenormalized theory 
of the arbitrary length scale introduced by renormal- 
ization the correlation functions 



of the VT order parameter satisfy the identity 



(20) 



(21) 



This identity translates with the help of the renormahza- 
tion scheme M) via the Gell-Mann-Low- Wilson functions 



(3iu) 
k{u) 
7(u) 



du 
dfj. 



91nr 



dliiZ 



(22) 



where the bare quantities are kept fixed while taking the 
derivatives, into the RGE 



du dr 2 ^ 



(23) 



Since the functions /3, k, 7 are entirely determined by the 
renormahzation factors, im particular by the Y^^\u), the 
RGE (|2^) is exactly equal to the RGE for the correlation 
functions of percolation theory. The RGE determines the 
scaling structure of a field theory. Thus, the VT has the 
same critical exponents as percolation. 



ly = [2 — k(u*)] 



77 = 7(u,) , 



(24) 



where is the infrared stable fixed point determined by 
/?(u*) = 0. The percolation exponents are known to third 
order e [|ll|. Introducing the correlation length ^ ^ \t\~^ 
and the order parameter exponent /? — v{d —2-1- ?7)/2, 
the asymptotic solution of the RGE ( |2^ ) can be written 
as 



(25) 



The RRN differs from percolation as long as the bond- 
resistance w is finite. The bond-resistance constitutes a 
further scaling variable and leads to an additional deriva- 
tive C^wd/dw, where C,{u) = /^^^j^|o, in the correspond- 
ing RGE 1^ . Defining the resistance crossover exponent 
by = i^[2 — C(u*)], the asymptotic solution for the cor- 
relation functions of RRN can be written as 

G(^^^^({x,0l,r,H 

= e-'^^/'^Ff ''''^({x/e, e/^/^i^}) . (26) 

Although the VT and percolation are showing the same 
scaling behavior, the universal scaling functions 

and limto^o f'^^'^^ may be different for both theories. 
This is due to the fact that these functions are not en- 
tirely determined by the superficially divergent parts of 
the corresponding vertex functions. The universal scaling 
functions cannot be calculated from the RGE. 



IV. ORDER PARAMETERS 

To explore the relation between VT and percolation 
further we now revisit the order parameters. We start 
with the order parameter for the RRN, Eg. (o ). Above 
the percolation point we deduce from Eq. (|2q ) and the 
constraint (0) that 



xp(^^^)(x)exp 



(x,( 



= Foo(t) / dx 



(27) 



Here Poo{t) ~ \t\^ is the percolation probability and 
i(;^'/'/'^p(RRN)^^^0/i.j.^) l-j^g probability distribution of 

the conductance Sqo = 1 / Roo from an arbitrary point on 
the percolating cluster to infinity. The digit 1 in the 
bracket denotes the limit of 1/n = (2A/)^^ for D 
0. In the n-state Potts model limit w 0, we retrieve 
from Eq. (|2^) the Potts order parameter M {9, r) = 
Poo(r)(<5,-o-l). 

The order parameter for the VT has an analogous 
form 0, |, I 



M(VT)(f,r) = (f](f))^^^=P^(r)|^' 



dx 



- 1 



(28) 



where = ('*''"'' ~ corresponds to the ra- 

dius of gyration of the replicas about a center of mass 
at f = ^ELii"^"^- ^^p^^'^H^^a) is the probability 
distribution of the inverse squares a of the localization 
lengths. Here, the digit 1 in the bracket denotes the limit 
of for n ^ 1. In general the two distributions 



P 



(VT) 



(x) and p 



(RRN) I 



x) are different but universal. A 



mean field calculation leads in both cases (for the RRN 
see Appendix 0) to Poo{t) = 2\T\/g and to the same 
integro-differential equation ||. || 




p{x) 



dx' p{x ~ x')p{x') , (29) 



with p{0) ~ p{oo) — 0. It is interesting that the Laplace 
transform p(z) = dx p{x) exp{— zx) and the differen- 
tial equation that follows from Eq. (ES), 



4zp"(z) = (l-25(z))p(z), 



(30) 



with p{0) ~ 1 and p(oo) — 0, was already introduced 
many years ago by Stephen and Stinchcombc fl^ in 
their mean field theory of RRN and determination of the 
conductivity of a Bethe lattice, respectively. 



V. CONCLUDING REMARKS 

We showed that the primitive divergences occurring in 
percolation and vulcanization are identical to arbitrary 
order in perturbation theory. Consequently both tran- 
sitions are governed by the same Gell-Mann-Low RGE. 
This RGE determines entirely the scaling structure of 
both transitions. Hence, the scaling behavior of physi- 
cally analogous quantities in both transitions is identi- 
cal. In particular, the VT is governed by the critical 
exponents of the percolation universality class. Quanti- 
ties which are not completely determined by superficial 
divergences cannot be calculated from the Gell-Mann- 
Low RGE. These quantities, e.g. scaling functions, may 
be different in both theories. 

Moreover, we compared the order parameter for the 
RRN and the VT. These two have an analogous form. 
The order parameter for the RRN involves a scaling func- 
tion which incorporates the distribution of the conduc- 
tance to infinity whereas the order parameter for the VT 
features the distribution of the inverse squares of local- 
ization length. In mean field approximation, the corre- 
sponding distribution functions turn out to be identical. 
To determine these distribution functions beyond mean 
field level is an interesting issue for future work. 

We would like to emphasize that the analysis of scaling 
properties presented in Sec. Ill is hardly feasible without 
employing the methods of renormalized local field the- 
ory. These methods provide a clear cut discrimination 
between scaling properties and other universal quantities 
like scaling functions. Furthermore, these methods allow 
to restrict attention to superficially divergent diagrams 
which simplifies the analysis tremendously. 

The intimate relationship between vulcanization and 
percolation seems plausible because macroscopic connec- 
tivity is the central issue in both systems. There are, 
on the other hand, striking distinctions between vulcan- 
ization and percolation. For example, common percola- 
tion models like the RRN live on some underlying lat- 
tice whereas vulcanization involves no lattice. Vulcan- 
ization as considered in this paper features an excluded 
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FIG. 1: Reduction of the one-loop self-energy diagram. 



volume interaction which is extraneous to percolation. A 
similarity between the excluded volume interaction and 
the lattice is, however, that both prevent overlap of the 
network constituents, monomers and bonds, respectively. 
Another obvious distinction is that the usual percolation 
involves a single ensemble, viz. the ensemble of the di- 
luted lattice configurations, whereas vulcanization fea- 
tures two ensembles: one pertaining to the thermal de- 
grees of freedom and one to the crosslink distribution. It 
turns out, though, that fluctuations of the crosslink dis- 
tribution play a more important role than thermal fluc- 
tuations do, at least as far as the connectivity aspects of 
the VT are concerned. 
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APPENDIX A: REDUCTION OF 
SUPERFICIALLY DIVERGENT DIAGRAMS 

In this appendix we demonstrate the simplicity of de- 
termining the one-loop contributions to the renormaliza- 
tion of the VT model. Consider the self-energy diagram 
constructed of the HRS (bold) propagators. Fig. 1. After 
decomposition in amorphous (thin) and liquid (dashed) 
propagators, we get the second line in Fig. 1, where sum- 
mation over the colors a, (3 is implied. We have shown 
that diagrams with different colors do not contribute to 
the primitive divergencies. Hence, we can set a = f3 in 
the third diagram. However, then the third diagram can 
be discarded because the wave vectors of the two propa- 
gators of the same color can not add to a HRS external 
wave vector at the vertices. Now we can replace all the 
propagators by the simple one (r-|-q^)~^ which results in 
the simple diagram appearing in the third line of Fig. 1. 
Therefore we get in dimensional regularization, after in- 
cluding the combinatorial factor 1/2, the following one- 
loop contribution p^~i°°P ^j^g two-leg vertex function 
r2 (note that the vertex functions Fjv are defined as the 
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A 




-3 




a 



= -2 




FIG. 2: Reduction of the one-loop vertex diagram. 

negative of the corresponding diagrams): 
,2 r 1 



r 



1 — loop 



T+ (q + k)2 

k2 



l-e/2 



6 



where e = 6 - d and = (An)'^^^ /T{1 + e/2). Renor- 
maUzation, Eq.(|), with Ta? ^ f^v = Z^^/^Tn and the 
identification = K'^ leads to the 1-loop result 



Zr 



Z - 



6(1 -s/2) 
^ f 2/ Y^^ 



e/2 



The e-poles are eliminated by choosing minimally 

Z=\ + ^ + 0(u2) , Zr = l + - + O(u') . (A3) 
be £ 

Now we consider the first order correction pi~'°°P 
to the vertex function T^. The three-leg diagram con- 
structed of the HRS (bold) propagators, Fig. 2, de- 
composes in amorphous and liquid (dashed) propagators 
shown in the second line of Fig. 2. Once more the sum- 
mation over the colors a, (3, and 7 is implied. We can set 
a = (3 — J in the diagrams to extract the primitive diver- 
gencies. However, then the third and fourth diagram can 
be discarded because the wave vectors of the propagators 
of the same color cannot add to a HRS external wave vec- 
tor at all the vertices. We are left the third line of Fig. 2 
where the propagators are replaced by (r + qf )^^. Thus, 
we find 



1 



r3({0}) = 2g3 

/q [T + q^ 

After renormalization we get 



2.g3T-^/2 



A, 



(A4) 



r^-'°°p({o},T) 



e/2 



(A5) 



It follows the remaining renormalization factor 



Zu = 1 



4U ^/ 9n 



(A6) 



All the renormalization factors, Eqs.(A^,A6) are known 
since a long time as the factors that rcnormalize the 
Potts model in the one state limit Denoting the 

logarithmic derivatives of the renormalization factors by 
7... = /i9 In Z.../3/i|Q, one derives easily the Gell-Mann- 
Low- Wilson functions of the RGE as 



f3{u) = (^-£ + 37-7„^i 



= ( - e+^u + Oiu^)]u, (A7) 



7(m) 



6 

7-7- 



+0("'): 



(A8) 
(A9) 



The fixed point as the asymptotic solution of the flow 
equation ldu{l)/dl — (3{u) is found to be = 2£/7 + 
0{e^). The critical exponents follow as = = 
-e/2l + 0{e'^) and v-^ = 2-k^(u*) = 2~5£/21 + 0(£2). 
These are the known percolation exponents. 



APPENDIX B: THE ORDER PARAMETER OF 
THE RRN 

In this appendix we discuss some properties of the or- 
der parameter M(0, t, w) of the random resistor network 
above the percolation point. We begin with a mean field 
consideration based on the saddle point equation steam- 
ing from Eq. (|^): 



x,( 



i«V^U(x,0) 



(Bl) 



The last term in this equation follows from the condi- 
tion that the variation must be done subject to the con- 
straint (||). We seek a spatially homogeneous solution, 
rotational invariant about 9 — Q 'm replica space. Thus, 
we make the ansatz (^(x, 9) = M{6) = m{f{9^)~l). Here 
denotes / a localized function with Jg, f{6''^) = jg,^ — n. 
From Eq. (Bl) we get, assuming to 7^ for r < 0, 



= [T-wVi)j{e')-T 



gm 



^ ^fie'r - 2f{9^) + 2 - - f{9'rj (B2) 

We separate this equation in its localized and its delocal- 
ized parts: 



■q2 



^r~wWl)f{e')-^-^(f{9'r-2f{e')) = 0(B3) 



0(B4) 
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Now we write the function / as a Laplace integral 



f{e')= / dtp{t)eM-m 

Jo 



(B5) 



Then we have for the term containing the replica space 
derivative 

= y dtp{t)(4:t^(P -2tDjcxp{-te^) 

^j^dtAeM-te')j^[t^p(t)), (B6) 

where the last equality holds in the limit D ^ 0, i.e., 
n — > 1. Furthermore we can deduce in this limit from 
the normalization of / that /p dtp{t) — 1 which leads 

in the limit n -> 1 toJsf{0Y = Ugfi(^^)? = 1- Thus, 
we obtain from Eq. ( |B4| ) that the mean field percolation 
order parameter is given for r < by m = — 2r/g = 
2 |t| /g. Using this result and equating the coefficients of 
exp(— t^) we get finally from Eq. (B3) 



(t'pit)) = Pit) - j dt' p{t')p{t - t') , (B7) 



Aw d f , 



which constitutes, after a rescaling, the integro- 
differential equation (|2^). A very good approxima- 
tive solution of this equation is given by p{x) — 



'exp(— l/4x) for X <C 1 and p{x) 



24(6a; 



3/5)exp(— 6x) for a; ^ 1 with a — 0.56925 and b — 
13.424 These asymptotic forms yield also very good 
approximations in the overlapping region a; ~ 1. 

The previous mean field consideration is valid in an ex- 
act sense only for dimensions d > 6. To get information 
about the behavior of the distribution p{x) on the spatial 
dimension d below 6, we calculate the mean macroscopic 
conductance to infinity Sqo- We consider spatial length 



scales large in comparison to the correlation length ^. In 
this regime the RRN above the percolation point, r < 0, 
can be considered as a homogeneous material of conduc- 
tivity ct{t). Instead of addressing Eqo directly we de- 
termine the resistance Rl of the RRN between in inner 
sphere of radius a 3> ^ and an outer one of radius L 3> a. 
The solution to this simple problem may be gleaned from 
many textbooks on electricity and magnetism. One finds 
that the voltage V behaves as a function of the current 
/ as 



V 



(j{d-2)Sd \ a 



1 



1 



L 



d-2 



RlI, (B8) 



where Sd is the surface of the unit sphere in d dimensions. 
Soo can now be obtained by taking the limit L ^ oo. The 
leading terms in this limit are 



S r — R, 



crSd 



{d - 2)a'^-2 fQj, d > 2 , 
(ln(i/a)) ford = 2,(B9) 

(2 - for d < 2 . 

Hence, the macroscopic conductance to infinity is fi- 
nite in the case of d > 2 and vanishes for d < 2. We 
conclude that the distribution of the conductance to in- 
finity p{x) must develop a 5-peak at a; = if d < 2 and 
the order parameter of the RRN vanishes. If we make a 
scaling ansatz 



Ei(r) = |r|*F(L/C,a/0 



(BIO) 



and compare with Eq.(B9), we find the well known ex- 



ponent t for the macroscopic conductivity (t(t) — \t\ 

t={d-2)v + (j) (Bll) 

in all cases. 
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